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Abstract 

We explore a class of effective field theory models of cosmic acceleration involving a metric and 
a single scalar field. These models can be obtained by starting with a set of ultralight pseudo- 
Nambu-Goldstone bosons whose couplings to matter satisfy the weak equivalence principle, as- 
suming that one boson is lighter than all the others, and integrating out the heavier fields. The 
result is a quintessence model with matter coupling, together with a series of correction terms 
in the action in a covariant derivative expansion, with specific scalings for the coefficients. After 
eliminating higher derivative terms and exploiting the field redefinition freedom, we show that 
the resulting theory contains nine independent free functions of the scalar field when truncated 
at four derivatives. This is in contrast to the four free functions found in similar theories of 
single- field inflation, where matter is not present. We discuss several different representations of 
the theory that can be obtained using the field redefinition freedom. For perturbations to the 
quintessence field today on subhorizon lengthscales larger than the Compton wavelength of the 
heavy fields, the theory is weakly coupled and natural in the sense of t'Hooft. The theory admits 
a regime where the perturbations become modestly nonlinear, but very strong nonlinearities lie 
outside its domain of validity. 
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1 Introduction and Summary 



1.1 Background and Motivation 

The recent discovery of the accelerating expansion of the Universe [U [2] has prompted many 
theoretical speculations about the underlying mechanism. The most likely mechanism is a cos- 
mological constant, which is the simplest model and is in good agreement with observational data 
[3] . More complicated models involve new dynamical sources of gravity that act as dark energy, 
and/or modifications to general relativity on large scales. A plethora of models have been postu- 
lated and explored in recent years, including Quintessence, K-essence [HE], Ghost Condensates 
[6J, DGP gravity [7], and f(R) gravity, to name but a few. See Refs. [ffl EQ E21 [IS E] for 
detailed reviews of these and other models. 

A common feature of the majority of dark energy and modified gravity models is that in the 
low energy limit, they are equivalent to general relativity coupled to one or more scalar fields, 
often called quintessence fields. Therefore it is useful to try to construct very general low energy 
effective quantum field theories of general relativity coupled to light scalar fields, in order to 
encompass broad classes of dark energy models. Considering dark energy models as quantum 
field theories is useful, even though the dynamics of dark energy is likely in a classical regime, 
because it facilitates discriminating against theories which are theoretically inconsistent or require 
fine tuning. 

A similar situation occurs in the study of models of inflation, where it is useful to construct 
generic theories using effective field theory. Cheung et al. [15] constructed a general effective field 
theory for gravity and a single inflaton field, for perturbations about a background Friedman- 
Robertson- Walker cosmology in unitary gauge. This work was later generalized in multiple 
directions |16t [T7] and has been very useful. An alternative approach to single field inflationary 
models was taken by Weinberg [18], who constructed an effective field theory to describe both 
the background cosmology and the perturbations. This theory consisted at leading order of a 
standard single field inflationary model with a potential, together with higher order terms in a 
covariant derivative expansion up to four derivatives. More detailed discussions of this type of 
effective field theory were given by Burgess, Lee and Trott |19j . 

When one turns from inflationary effective field theories to quintessence effective field theories, 
the essential physics is very similar, but there are three important differences that arise: 

• First, the hierarchy of scales is vastly more extreme in quintessence models. The Hubble 
parameter H is typically several orders of magnitude below the Planck scale m p ~ 10 28 eV 
in inflationary models, whereas for quintessence models Hq ~ 10~ 33 eV is ~ 60 orders of 
magnitude below the Planck scale. Quintessence fields must have a mass that is smaller 
than or on the order of Hq. It is a well-known, generic challenge for quintessence models 
to ensure that loop effects do not give rise to a mass much larger than Hq. Because of 
the disparity of scales, this issue is more extreme for quintessence models than inflationary 
models. 

• In most inflationary models, it is assumed that the dynamics of the Universe are dominated 
by gravity and the scalar field (at least until reheating). By contrast, for quintessence 
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models in the regime of low redshifts relevant to observations, we know that cold dark 
matter gives an 0(1) contribution to the energy density. Therefore there are additional 
possible couplings and terms that must be included in an effective field theory. 

For any effective field theory, it is possible to pass outside the domain of validity of the theory 
even at energies E low compared to the theory's cutoff A, if the mode occupation numbers N 



are sufficiently large (see Sec. 5.2 below for more details). This corresponds to a breakdown 
of the classical derivative expansion. For quintessence theories, mode occupation numbers 
today can be as large as N ~ (m p /Ho) 2 and it is possible to pass outside the domain of 
validity of the theory. By contrast in inflationary models, this is less likely to occur since 
mode occupation numbers for the perturbations are not large before modes exit the horizon. 
Thus, the effective field theory framework is less all-encompassing for quintessence models 
than for inflation models. This issue seems not to have been appreciated in the literature 
and we discuss it in Sec. 15.21 below. 

Several studies have been made of generic effective field theories of dark energy. Crem- 
inelli, D'Amico, Noreha and Vernizzi [20] constructed a the general effective theory of single-field 
quintessence for perturbations about an arbitrary FRW background, paralleling the similar con- 
struction for inflation [15]. Park, Watson and Zurek constructed an effective theory for describing 
both the background cosmology and the perturbations, following the approach of Weinberg j!8j 
but generalizing it to include couplings to matter [21]. 

The two approaches to effective field theories of quintessence - specialization to perturba- 
tions about a specific background, and maintaining covariance and the ability to describe the 
dynamics of a variety of backgrounds - are complementary to one another. The dynamics of the 
cosmological background FRW solution can be addressed in the covariant approach of Weinberg, 
but not in the background specific approach of Creminelli et al., which restricts attention to the 
dynamics of perturbations about a given, fixed background. On the other hand, a background 
specific approach can describe a larger set of dynamical theories for the perturbations than can 
a covariant derivative expansion^ 

1.2 Approach and Assumptions 

The purpose of this paper is to revisit, generalize and correct slightly the covariant effective field 
theory analysis of Park, Watson and Zurek [21]. Following Weinberg and Park et al., we restrict 
attention to theories where the only dynamical degrees of freedom are a graviton and a single 
scalar. We allow couplings to an arbitrary matter sector, but we assume the validity of the 
weak equivalence principle, motivated by the strong experimental evidence for this principle. We 
assume that the theory consists of a standard quintessence theory coupled to matter at leading 

1 To see this, consider for example a term in the Lagrangian of the form f((f>)(X7(f>) 2n , where <f> is the quintessence 
field. Such a term would be omitted in the covariant derivative expansion for sufficiently large n. However, upon 
expanding this term using (f> — (fio + 8<j), where 4>o is the background solution, one finds terms ~ (V<?!>o) 2n ~ 2 (V<5<?!>) 2 
which are included in the Creminelli approach of applying standard effective field theory methods to the pertur- 
bations. 
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order in a derivative expansion, with an action of the form 

S[g a pAM = J d A x^~g | - l -(Vcft) 2 -U(<f>) j + S m [e^^Vm] . (1.1) 

Here ip m denotes a set of matter fields, and m p is the Planck mass. The factor e a ^ in the matter 
action provides a leading-order non- minimal coupling of the quintessence field to matter, in a 
manner similar to Brans-Dicke models in the Einstein frame. 
Our analysis then consists of a series of steps: 

1. We add to the action all possible terms involving the scalar field and metric, in a covariant 
derivative expansion up to four derivatives. We truncate the expansion at four derivatives, 



as this is sufficient to yield the leading corrections to the action (1.1). As described by 



Weinberg |18| there are ten possible terms, with coefficients that can be arbitrary functions 



of <ft [see Eq. (2.3) below]. Section 5.1 below describes one possible justification of this 
covariant derivative expansion from an effective field theory viewpoint, starting from a set 
of ultralight pseudo Nambu-Goldstone bosons (PNGBs). It is likely that the same expansion 
can be obtained from other, more general starting points. 

2. We allow for corrections to the coupling to matter by adding to the metric that appears in 
the matter action all possible terms involving the metric and (ft allowed by the derivative 



expansion, that is, up to two derivatives. There are six such terms [see Eq. (2.4) below.] 
We also add to the action terms involving the stress energy tensor T^ u of the matter fields, 



up to the order allowed by the derivative expansion using T^ u ~ miG^ [see Eq. (2.3) 



below]. Including such terms in the action seems poorly motivated, since a priori there is 
no reason to expect that the resulting theory would respect the weak equivalence principle 
(see Appendix [B]) . However we show in Appendix |B| that the weak equivalence principle is 
actually satisfied, to the order we are working to in the derivative expansion. In addition, 
all the terms in the action involving T^ v can be shown to have equivalent representations 
not involving the stress energy tensor, using field redefinitions (see Appendix [B|) . 

3. The various correction terms are not all independent because of the freedom to perform field 
redefinitions involving (ft, g^v and the matter fields, again in a derivative expansion. In Sec. 
[3]we explore the space of such field redefinitions, finding eleven independent transformations 
and tabulating their effects on the coefficients in the action (see Table [T] below) . 

4. Several of the correction terms that are obtained from the derivative expansion are "higher 
derivative" terms, by which we mean that they give contributions to the equations of 
motion which involve third-order or higher-order time derivatives of the fields^] Normally, 

2 The precise definition of higher derivative that we use, which is covariant, is that an equation will be said not 
to contain any higher derivative terms if there exists a choice of foliation of spacetime for which any third-order or 
higher-order derivatives contain at most two time derivatives. Theories which are higher derivative in this sense 
are generically associated with instabilities (Ostragradski's theorem) [22 , although the instabilities can be evaded 
in special cases, for example f(R) gravity. For most of this paper (except for the Chern-Simons term), a simpler 
definition of higher derivative would be sufficient: a term in the action is "higher derivative" if it gives rise to terms 
in the equation of motion that involve any third-order or higher order derivatives. 
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such higher derivative terms give rise to additional degrees of freedom. However, if they 
are treated perturbatively (consistent with our derivative expansion) additional degrees of 
freedom do not arise. Specifically, one can perform a reduction of order procedure on the 
equations of motion [23j [231 [22] j substituting the zeroth-order equations of motion into the 
higher derivative terms in the equations of motion to eliminate the higher derivatives^} We 
actually use a slightly different but equivalent procedure of eliminating the higher derivative 
terms directly in the action using field redefinitions^] (see Appendix [c| . 

Weinberg [TH] and Park et al. [H] use a slightly different method, consisting of substituting 
the leading order equations of motion directly into the higher derivative terms in the action. 
This method is not generally valid, but it is valid up to field redefinitions that do not involve 
higher derivatives, and so it suffices for the purpose of attempting to classify general theories 
of dark energy (see Appendix |C]) . 

5. Another issue that arises with respect to the higher derivative terms is the following. Is 
it really necessary to include such terms in an action when trying to write down the most 
general theory of gravity and a scalar field, in a derivative expansion? Weinberg |18) 
suggested that perhaps a more general class of theories is generated by including these 
terms and performing a reduction of order procedure on them, rather than by omitting 
them. However, since it is ultimately possible to obtain a theory that is perturbatively 
equivalent to the higher derivative theory, and which has second order equations of motion, 
it should be possible just to write down the action for this reduced theory. In other words, an 
equivalent class of theories should be obtained simply by omitting all the higher derivative 
terms from the start. We show explicitly in Sec. [4] that this is the case for the class of 
theories considered here. 

6. We fix the remaining field redefinition freedom by choosing a "gauge" in field space, thus 



fixing the action uniquely (see Sec. 4.2). 
1.3 Results and Implications 



Our final action is [Eq. (4.5) below] 
S = f d 4 x^g I - t 



l 

E R - _(V0) 2 - U{4>) \ + S m [^g a p, 
+ e J d 4 x^|ai(V</») 4 + 6 2 T(V0) 2 + ciG^V^V^ 

+ d 3 {R 2 - 4R^R^ + R^pR^o) + die^C^Cxpap + exT^T^ + e 2 T 2 + . . . j . 



(1.2) 



3 This is more general than requiring the solutions of the equation of motion to be analytic in the expansion 
parameter, as advocated by Simon [26) : see Ref. [25] . 

4 This procedure is counterintuitive since normally field redefinitions do not change the physical content of a 
theory; here however they do because the field redefinitions themselves involve higher derivatives. 
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Here the coefficients a±, &2 etc. of the next-to-leading order terms in the derivative expansion 
are arbitrary functions of 0, and the ellipsis . . . refers to higher order terms with more than four 
derivatives. The corresponding equations of motion do not contain any higher derivative terms. 
This result generalizes that of Weinberg [18] to include couplings to matter. 
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Figure 1: The parameter space of fractional density perturbation Sp/ p for perturbations to the quintessence 



Geld, and cutoff scale M for the effective held theory, illustrating the constraint [1.3) on the domain of 



validity. Near the boundary of the domain of validity the higher derivative terms in the action are 
potentially observable, this is labeled the "interesting regime". Further away from the boundary the 
higher derivative terms are negligible and the theory reduces to a standard quintessence model with a 
matter coupling. 



We can summarize our key results as follows: 

• The most general action contains nine free functions of 4>: U, a, oi, £>2> ci, cfe, d^, e\, e2, as 
compared to the four functions that are needed when matter is not present |18| . 



There are a variety of different forms of the final theory that can be obtained using field 
redefinitions. In particular some of the matter-coupling terms in the action can be re- 
expressed as terms that involve only the quintessence field and metric. Specifically, the 
term T(V(/>) 2 term could be eliminated in favor of □0(V^) 2 , the (V0) 4 could be eliminated 
in favor of a term T^Vu^V^, or the G^Vu0V y term could be eliminated in favor of a 



term T^V v <t> (see Sec. 4.2). 



As mentioned above, one obtains the correct final action if one excludes throughout the 
calculation all higher derivative terms. 

The final theory does contain terms involving the matter stress-energy tensor. Nevertheless, 
the weak equivalence principle is still satisfied (see Appendix [B]) . It is possible to eliminate 
the stress-energy terms, but only if we allow higher derivative terms in the action (where it 
is assumed that the reduction of order procedure will be applied to these higher derivative 
terms). Thus, for a fully general theory, one must have either stress-energy terms or higher 



derivative terms; one cannot eliminate both (see Sec. 4.2). 
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We can estimate how all the coefficients a\ etc. scale with respect to a cutoff scale M for 



an effective field theory as follows (see Sec. 5.1). We assume that several ultralight scalar 



fields of mass ~ Hq arise as pseudo Nambu-Goldstone bosons (PNGBs) from some high 
energy theory \17\ I28j . and are described by a nonlinear sigma model at low energies. We 
then suppose that all but one of the these PNGB fields have masses M that are somewhat 
larger than ~ Hq, and integrate them out. This will give rise to a theory of the form 
discussed above for the single light scalar, where the higher derivative terms are suppressed 
by powers of M. The scalings for each of the coefficients in the action are summarized in 
Table [3| We find that the fractional corrections to the cosmological dynamics due to the 
higher derivative terms scale as Hq/M 2 , as one would expect. 

Finally, we can use these scalings to estimate the domain of validity of the effective field 



theory (see Sec. 5.2). We find that cosmological perturbations with a density perturbation 



5p in the quintessence field must have a fractional density perturbation that satisfies 



« T72 • (1-3) 



5p M 2 

n^~R 2 
P H 

Thus perturbations can become nonlinear, but only modestly so, if M is close to Hq. The 
parameter space of fractional density perturbation 5p/p and cutoff scale M is illustrated in 
Fig. [Tj In addit ion there is the standard constraint for derivative expansions 

E<^M (1.4) 

where E~ x is the length-scale or time-scale for some process. We show in Fig. [2] the two 



constraints (1.3) and (1.4) on the two dimensional parameter space of energy E and mode 



occupation number N. 

Finally, in Appendix [D] we compare our analysis to that of Park, Watson and Zurek 
who perform a similar computation but in the Jordan frame rather than the Einstein frame (see 
also Ref. |29j). The main difference between our analysis and theirs is that they use a different 
method to estimate the scalings of the coefficients, and as a result their final action differs from 
ours, being parameterized by three free functions rather than nine. 

2 Class of Theories Involving Gravity and a Scalar Field 

As discussed in the Introduction, our starting point is an action for a standard quintessence model 
with an arbitrary matter coupling, together with a perturbative correction which consists of a 
general derivative expansion up to four derivatives. The action is a functional of the Einstein- 
frame metric g a a, the quintessence field (ft, and some matter fields which we denote collectively 
by ip m : 

S[5 Q/3 ,0, Vm] = So[g a p,4>] + e5i[g a( g,0,r a ^(^ m )] + S m \g a p, ipm] + 0(e 2 ). (2.1) 
Here S m is the action for the matter fields, and the quantity e is a formal expansion parameter. 



We will see in Sec. 5.1 below that e can be identified as proportional to M , where M is a 
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Figure 2: The domain of validity of the effective field theory in the two dimensional parameter space 
of energy E per quantum of a mode of the quintessence Geld, and mode occupation number N. The 
cutoff scale M must he larger than the Hubble parameter Hq in order that the background cosmology lie 
within the domain of validity. Perturbation modes on length-scales that are small compared to Hq 1 but 
large compared to M _1 can be described, but only if the mode occupation number and fractional density 



perturbation are sufficiently small. See Sec. 5.2 for details. 



cutoff scafe or the mass of the lightest of the fields that have been integrated out to obtain the 
low energy action. Equivalently, e counts the number of derivatives in our derivative expansion, 
with e n corresponding to 2n derivatives. The notation in the second term indicates that the 
perturbative correction S\ to the action can depend on the matter fields, but only through their 
stress energy tensor T a p (defined in Appendix [A|. Explicitly we have 

,2 



So 



d 4 Xy 



-fj 



mz 1 
_^__( V 0) 2 -[/(0) 



(2.2) 



and PSIEU 



Si 



- 5 ja^) 4 + a 2 D0(V</>) 2 + a 3 (n^) 2 

+ b 6 RT + b 7 T + ciG^V^V^ + c 2 R(V(f>) 2 + c^RUcp 



+ dxR 2 + d 2 R^R 



d 3 (R 2 - AR^R^ + R^apR^P) 



+ d^C^Cx^ + exT^T^ + e 2 T' 



(2.3) 
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Here C^xp is the Weyl tensor and e^ u ^ p is the antisymmetric tensor (our conventions for these 
are given in Appendix |A| . There are additional terms with four derivatives that one can write 
down, but all such terms can be eliminated by integration by parts. Finally, the metric which 
appears in the matter action S m in Eq. (2.1) is given bjj^] 



9v» =e a g^ + ee a [faV „<f> + /3 2 (Vc/>) 2 ^ + faUfa^ + /^V^ + faR^ + faRg^] 

+ 0(e 2 ). (2.4) 

All of the coefficients o«, bi, Cj, di, ej, fa and a are arbitrary functions of <f>. 

Let us briefly discuss each of the perturbative terms. The terms with coefficients a% are 
corrections to the kinetic term of the scalar field. The bi and fa terms are couplings between 
the scalar field and the stress-energy tensor, or between curvature and the stress-energy ten- 
sor. The Ci terms are kinetic couplings between the scalar field and gravity. The di terms are 
quadratic curvature terms, which we have chosen to write as an R 2 term, an Rp V R^ u term, and 
the Gauss-Bonnet term. Any constant piece of the coefficient d% is a topological term and may 
be omitted. The term d^ is the gravitational Chern-Simons term, which may be excluded if one 
wishes to introduce parity as a symmetry of the theory, and again, any constant component of 
c?4 is topological and may be omitted. Finally, the terms are quadratic in the stress-energy 
tensor. 



Note that several of the terms in the action (2.3) are "higher derivative" terms, that is, they 



give rise to contributions to the equations of motion containing derivatives of order three or 
higher. The specific terms are those parameterized by the coefficients a^, 63, . . . , &g, c 2 , C3, d\, d 2 
and fa,...,fa. As discussed in the Introduction and in Appendix [Cj we will choose to define our 
theory by treating these terms perturbatively, which excludes the extra degrees of freedom and 
instabilities that are normally associated with higher derivative terms. 



We also note that the theory (2.1 ) satisfies the weak equivalence principle, to linear order in e, 
as we show in Appendix|B] That is, objects with negligible self-gravity with different compositions 
all experience the same acceleration. It is not a priori obvious that the principle should be satisfied 
since, as we show in Appendix [E] violations of the principle generically arise whenever the matter 



stress energy tensor appears explicitly in the gravitational action, as in Eq. (2.1) 

3 Transformation Properties of the Action 



The description of the theory provided by Eqs. (2.1 ) - (2.4) is very redundant, in part because of 



the freedom to perform field redefinitions. In this section we derive how the various coefficients 



in the action (2.1) are modified under various transformations. In the next section we will use 
these transformation laws to derive a canonical representation of the theory, involving only nine 
free functions. 



5 We call this metric the Jordan frame metric, in an extension of the usual terminology which applies to the case 



when the relation (2.4| between the two metrics is just a conformal transformation. 
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3.1 Expansion of the Matter Action 



Consider first the perturbative terms parameterized by f3±, . . . in the definition (2.4) of the 



Jordan metric g a p, which appears in the matter action S m [g a p,ip m ]. Using the definition (A.l) 



of the stress-energy tensor, we can eliminate these terms in favor of terms in the action involving 



T a R. Specifically we have from Eq. (A.l) that 



S m [e a (g^ + 5g^ v )^ m ] = S m [e a g^ ^m] + 2 J d 4 x^e 2a T^5 9tlu + 0(6g< 

Choosing 

Sg^ = e[faV^V u (/) + fa{V4>fg^ v + fa^toi* + AV^V^ + faR^ v + faRg^] 



(3.1) 



(3.2) 



then gives a transformation of the action (2.1) characterized by the following changes in the 
coefficients: 



(3.3) 



5pi = 


-fa, 


Sbi = 


h 2a fa 


6fo = 


-fa, 


5b 2 = 


h 2a fa 


<% = 


-fa, 


5b 3 = 


h 2a fa 




-fa, 


5b 4 = 


h 2a fa 


<5/3 5 = 


-fa, 


5h = 


\e 2a fa 


8fo = 


-fa, 


5b 6 = 


h 2a fa 



Here the parameters can be arbitrary functions of <j>. Similarly choosing 5g^ u = eag^ u gives a 
transformation characterized by 



5a 



-eat 



5b 7 



l„2a 



e a. 



(3.4) 



3.2 Field Redefinitions Involving just the Scalar Field 

Consider a perturbative field redefinition of the form 

4> = V> + £7, (3.5) 
where the quantity 7 can in general depend on any of the fields and their derivatives. To leading 



order in e, the change in the action (2.1) is then proportional to the zeroth-order equation of 

(3.6) 



motion (5.10b) for cf>. Relabeling ifi as (j), the change induced in the action is 

1 



5S 



d x^/^g^ 



U' + ^e 2a a'T 



There are three special cases that will be useful: 
1. First, choose 



tp + ea\T, 



(3.7) 
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where o\ is an arbitrary functional of ijj, and T is the trace of the stress-energy tensor. 



Substituting this into Eq. ( |3.6[ ) and comparing with the general action (2.3), we find the 
following transformation law for the coefficients: 

5bs = o"i, <5&7 = —U'ai, 

5e 2 



2. Second, we use the field redefinition 

<P = ^ + ea 2 [U^ + U'{^)\. (3.9) 

Here the second term in the square bracket is included in order to maintain canonical 
normalization of the scalar field, that is, to avoid generating terms in the action of the form 
f(4>) (V</>) 2 . The resulting transformation law is 



5a 3 = a 2 , 5b 3 = \e 2a a'a 2 . 

5b 7 = \a'e 2a U'a 2 , 6U = e(U') 2 a 2 . 



2 (3.10) 



3. Third, consider the field redefinition 

«/, = V> + «T3-e4(VV) 2 /^ (3.11) 

where 03 is a function of tp and again the particular combination of terms is chosen to 
maintain canonical normalization. Substituting into Eq. (3.6), performing some integrations 
by parts and comparing with Eq. (2.3) gives the transformation law 

Sa 2 = -a' 3 /U', 6b 2 = -\e 2a a' a'JV , 

5b 7 = \e 2a a'a Zl 5U = eU'a 3 . [ ' ' 

Note that this transformation is not well defined in general in the limit U' — > 0, because 
of the factors of 1/U'. However, it is well defined in the limit U' — > 0, o" 3 — >• with cr' 3 /U' 
kept constant. 



3.3 Field Redefinitions Involving the Metric 

We now consider a more general class of field redefinitions, where in addition to redefining the 
scalar field via Eq. (3.5), we also perturbatively redefine the metric via 

g a p = 9a/3 + eF a 0. (3.13) 

Here the quantity F a p can depend on if), g a p, their derivatives and the stress energy tensor. The 
corresponding change in the action is proportional to the equation of motion (5.10a). Relabeling 
9a.fi as 9a/3 and ip as 0, the total change in the action is 



5S 



+ e d x^—gj 



mlG afi + V^V'V - \{V4>fg aP - Ug al3 + e 2a T c 



U<j)-U' + -e 2a a'T 



(3.14) 



'Because we are working to linear order in e, it does not matter whether we take o\ to be a function of </> or of 
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Coeff. 
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01 
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0-3 
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0"9 


<7lO 


011 


ai 
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□</>(V0) 2 
























03 t 
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T^V^V v <f> 
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T(V0) 2 
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64 t 


























^5 t 


























h t 


RT 








* 
















b 7 


T 








* 








* 








ci 


















* 








C2 t 


R(Vcj)) 2 
























C3 t 


RD(j> 
























di t 


R 2 
























d 2 t 


























d 3 


Gauss-Bonnet 
























di 


Chern-Simons 
























ei 


























e2 




























U (potential) 

























Table 1: This table shows which of the terms in our action (2.2) are affected by each of the eleven field 
redefinitions (3.7) - (3.29) that are parameterized by the functions <7i (</>), . . . , au((f>). The columns repre- 
sent the redefinitions, and the rows represent terms. Daggers f in first column indicate "higher derivative" 
terms, that is, terms that give contributions to the equations of motion containing time derivatives of 
higher than second order. Stars * indicate that the coefficient of that row's term is altered by that col- 
umn's field redefinition. We omit the coefficients a and /?i , . . . , /3g since those coefficients are degenerate 
with b x , . . . , b? by Eqs. JO) and pT ). 
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Note that this formula includes the effect of the change in the Jordan frame metric (2.4) caused 
by the transformation (3.13). We now consider seven different transformations of this type: 



4. The first case is a change to the metric proportional to Rg a p. In order to maintain canonical 
normalization of both the metric and the scalar field, that is, to avoid terms of the form 
f((f))(V(J)) 2 and f((f>)R, we need the following combination of terms in the field redefinition: 



mz 



Qa/3 =g a /3 ~ 2ea' A y-j^R + AJ g a[ s, 
4> =i> + 4ecr 4 , 



(3.15a) 
(3.15b) 



for some function a 4 (ip). Substituting into Eq. (3.14), performing some integrations by 
parts and comparing with Eq. (2.3) we obtain for the transformation law 

SC2 



5b 7 = 2e 2a a'a 4 - 4e 2 V 



U °4> 



1- 



5U = 4e [U'a 4 - AUa' 4 ] 



5ba 



U °4> 
e 2a 2 I 



(3.16) 



5. Next consider changes to the metric proportional to R a /3. In order to maintain canonical 
normalizations we use the following combination of terms in the field redefinition: 

m 2 

9a/3 =<?a/3(l - SecTg) - 2e-^-a' 5 R a/3 , (3.17a) 
(/>=V> + ecT5, (3.17b) 

for some function 05 (-0)- This gives the transformation law 

5b 7 = le 2a a'a 5 - e 2a aL 5a = -^V 



[T"5> 

5d x = -pa> 5 , 5d 2 = fa> 5 , (3.18) 

5b 5 = -^e 2a a' 5 , 5U = e[U'a 5 - AUa' 5 ] . 



6. The next case is a change to the metric proportional to (V0) 2 <? Qj g. To maintain canonical 
normalization of the scalar field, we need in addition a change to the scalar field, with the 
combined transformation being 

a' 

9a/3 =g a p ~ 2ejj(\7ip) 2 g af s, (3.19a) 
</>=^ + 4e<7 6 , (3.19b) 

for some function o§. The resulting transformation law for the coefficients is 

<5ai = a'jU, 5b 2 = -e 2a a' 6 /U, 

5b 7 = 2e 2a a'a 6 , 5c 2 = -a' & m 2 p /U, (3.20) 
SU = AeU'a 6 . 
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7. Next consider changes to the metric proportional to □^ Q , ( g. The required form of field 
redefinition that preserves canonical normalization of (p is 

9a/3 =9a/3 + 2ea 7 U^g af3 , (3.21a) 
(j) =i\> + 4eUa 7 , (3.21b) 

for some function a 7 . The coefficients in the action then change according to 

5a 2 = -a 7 , 5b 3 = e 2a a 7 , 

5b 7 = 2e 2a a'Ua 7 , 5c 3 = m 2 p a 7 , (3.22) 
SU = AeUU'a 7 . 

8. The fifth case is a change to the metric proportional to V a (j)V The required form of 
field redefinition that preserves canonical normalization of 4> is 

9a/3 =9a/3 - 2e jjV a tpV pip , (3.23a) 
(f)=ip + ecr 8 , (3.23b) 

for some function as- The coefficients in the action then change according to 

5ai = -a' 8 /(2U), 5h = -e 2a a' 8 /U, 

Sb 7 = le 2a a'a 8 , <5ci = m 2 p a' 8 /U, (3.24) 

SU = eU'a s . 

9. Next consider a change in the metric proportional to Vq.V/30. To preserve canonical nor- 
malization of (p we use the redefinitions 

g a /3 =9ai3 + 2ea 9 V a Vi3ip, (3.25a) 
<j)=^ + eUa 9 , (3.25b) 

for some function as- The coefficients in the action then change according to 

5ai = -7>a' 9 , 5a 2 = -cr 9 , 

Sb 4 = e 2a a 9 , 5b 7 = \e 2a a'Ua 9 , (3.26) 
5c\ = m 2 a' g , 5U = eUU'ag. 

10. A simple case is when the change in the metric is proportional to Tg a p, for which no change 
to the scalar field is required. The redefinition is 

9ap =9a/3 + 2ea 10 Tg a/3 , (3.27) 
for some function aio. The transformation law for the coefficients is 

5b 2 = —aio, 5b 7 = —AaioU, 



i< ;l ™2 r ( 3 - 28 ) 



5e 2 = e 2a a w , 5b 6 = mla w . 
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11. Similarly, no transformation to the scalar is required for the case of a change in the metric 
proportional to T a g. The redefinition is 

g a /3 =g a (3 + 2e<7nT a/3 , (3.29) 

for some function an, and the corresponding transformation law is 

Sbi = an, 5b 2 = -\an, 

Sb 7 = -a n U, 5e 1 = e 2a a lu (3.30) 
5b 5 = -m 2 p an, 5b 6 = \m 2 p an- 



The eleveiQ field redefinitions " (^29) 

are summarized in Table [TJ which shows which 



coefficients are modified by which transformations. 



4 Canonical Form of Action 



In this section, we derive our final, reduced action (1.2) from the starting action (2.1), using 
the transformation laws derived in Sec. [3l There is some freedom in which terms we choose 
to eliminate and which terms we choose to retain. We choose to eliminate all terms that give 
higher derivatives in the equations of motion, so that the final theory is not a "higher derivative" 
theory. However, even after this has been accomplished, there is still some freedom in how the 



final theory is represented. We discuss this further in Sec. 4.2 below. The order of operations in 
the derivation is important, since we need to take care that terms which we have already set to 
zero are not reintroduced by subsequent transformations. Table [2] summarizes our calculations 
and their effects on the coefficients in the action at each stage in the computation. 

4.1 Derivation 

The steps in the derivation are as follows: 



1. Elimination of Derivative Terms in the Jordan Frame Metric: The transformation (3.3) 



can be used to eliminate all of the terms involving derivatives in the Jordan frame metric 



(2.4), which are parameterized by the coefficients f3\, . . . ,f3§. This changes the coefficients of 
the terms in the action that depend linearly on the stress energy tensor, namely b\, . . . , be. 
As discussed in Appendix [E] these terms involving the stress-energy tensor look like they 
might violate the weak equivalence principle, but in fact they do not. 

Elimination of Higher Derivative, Quadratic in Curvature Terms: We next consider the 
terms in the action that are quadratic functions of curvature, whose coefficients are d\, 
d%, d^ and d±. The Chern-Simons term (d^) and the Gauss-Bonnet term (d^) give rise 
to well behaved equations of motion (in the sense that they not increase the number of 
degrees of freedom), so we do not attempt to eliminate these terms. By contrast, the terms 



7 We could also consider a twelfth redefinition given by g a p = g a p{)- — 2ea' 12 ), cj> — ip + e<7i2 — empCr^R/V . 
However this redefinition is not independent of the first eleven; the same effect can be achieved by choosing 
«n = -e 2a a' 12 /U', a :i = -a 12 , a 7 = a' 12 /U', a w = e 2a a'a' 12 /(2U'). 
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Table 2: This table shows the progression of manipulations we make in this section. The second column 



on the left lists the various terms in the action (2.3), or in the Jordan-frame metric (2.4). The first column 
lists the corresponding coefficients; daggers f indicate higher derivative terms. The numbers in the first 
row along the top refer to the numbered steps in the derivation in Sec. |4.1| The second row shows which 
transformation functions are used in each step. In the table, a star ★ indicates that the corresponding row's 
term receives a contribution from the corresponding column's reduction process, while — > indicates that 
the term has been eliminated. The check marks / in the last column indicate the remaining terms that 



are non-zero in the final action (4.5). Finally, circles o in the last column indicate terms that are nonzero 



in alternative forms of the final action obtained using the transformations (3.11) or (3.23), as discussed in 
Sec.|42l 
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proportional to the squares of the Ricci scalar and Ricci tensor, parameterized by d\ and 
d,2, do increase the number of degrees of freedom. We can eliminate these terms by using 
the transformations (3.15) and ( |3.17| ), with parameters chosen to be 



0"4 



#i^(di + d 2 /2), 

rat 



0-5 



U 



mz 



do 



(4.1) 



These transformations will then modify the coefficients 65, b^, 67, c\ and c 2 , as well as the 
potential U (see Table [I]). 

Elimination of some of the Linear Stress-Energy Terms: We next turn to terms which 
depend linearly on the stress-energy tensor, parameterized by 61,..., b§. First, we can 
eliminate the term b^T^V „V v 



by using the transformation (3.25) with 09 



This gives rise to changes in the coefficients ax, 02, 67, c\ as well as to the potential U. 
Second, we can eliminate the terms parameterized by 65 and b§ by using the transformations 
(3.27) and (3.29) with the parameters axo = — (&6 + bs/2)/rn 2 ,, o~xi 



bs/rrip. This changes 



the coefficients 61, 62, 67, ex and e 2 . 



4. Elimination of Kinetic Coupling of the Scalar to Curvature: We next focus on the terms 
which kinetically couple the scalar field to gravity, namely G^V ^(fiS/ u 4>, R(V(j)) 2 and R\3(p. 
The first of these does not produce higher derivative terms in the equation of motion, so 
we focus on the remaining two terms, which are parameterized by C2 and C3. These terms 



can be eliminated using the transformations (3.19) and (3.21), with the parameters chosen 
to be 

, U 



o~6 



mt 



rC2, 



07 



C3 



mt 



(4.2) 



These transformations then give rise to changes in the coefficients eti, 02, 62; ^3) &7 as well 
as to the potential U. 

5. Elimination of some of the Corrections to Scalar Field Kinetic Term: Our action includes 
three corrections to the scalar kinetic term, parameterized by ax, 02 and 03. Of these, only 
term 03 contributes higher order derivatives to the equations of motion. We eliminate this 



term, and also the term 02, by using the transformations (3.9) and (3.11) with 



a 2 - -</;;• rr :] - j a<j)U' 'a 2 . (4.3) 
This gives rise to corrections to the coefficients b 2 , 63 and 67 and to the potential U. 
6. Elimination of some Kinetic Couplings of the Scalar to Stress-Energy: We next turn to 



the term biT^Vu^Vj/^. We can eliminate this using the transformation (3.23) with the 
parameter choice 



or 8 = J d^bxUe- 2a . (4.4) 
This gives rise to changes in the coefficients ax, by, c\ and U, from Table [T] We can also 



eliminate the term b^TOcj) by using the transformation (3.7) with o"i = —63. This changes 
the coefficients e 2 and 67. 
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7. Elimination of Trace of Stress-Energy Tensor Term: The last step is to re-express the term 
67T in terms of an 0(e) correction to the conformal factor e a by using the transformation 



(3.4) with a 



-2e~ 2a b 7 . 



4.2 Canonical Form of Action and Discussion 



Applying the parameter specializations derived above to the action (2.1) we arrive at our final 
result: 



S 



d £1/— g 



1 

l R _ _ (V 0) 2 - U{4>) \ + S m [e a ^g a(3 ^ m } 
+ e j dtxy/^gLn (V0) 4 + 6 2 T(Vc/>) 2 + ciG^V^V^ 

+ d 3 [R 2 - AR^R^ + R^pR^f) + d A e^C^C XpaP + e^T^ + e 2 T 2 }. 



(4.5) 



This action contains nine free functions of <p: U, a, a±, 6 2 , ci, d%, d±, ei, e 2 . The corresponding 
equations of motion do not contain any higher derivative terms and are presented in Appendix 

m 



Our final result (4.5) can be re-expressed in a number of equivalent forms: 



First, the term 6 2 T(V(/)) 2 in the action can be eliminated in favor of a term proportional 
to e 2a /3 2 (V<^) 2 5f /JI / in the Jordan frame metric (2.4) using the transformation (3.3). As 



discussed in Appendix [B] the latter representation makes explicit that the weak equivalence 
principle is satisfied. 

The term 6 2 T(V</>) 2 could also be eliminated in favor of a term a 2 D0(V(/>) 2 , using the 
transformation (3.11) parameterized by o 3 , as long as a' 7^ (j^J The dynamics of a scalar 
quintessence field with kinetic terms of the latter type have recently been explored in detail 
in Ref. [30], who called the mixing of the scalar and metric kinetic terms in the equations 
of motion "kinetic braiding". The representation of this term as a 2 n</>(Vi^) 2 has some 
advantages for cosmological analyses: in this representation the dynamics of the term are 
confined to the scalar sector, while in the 6 2 representation they are coupled to matter. 

The term ai(V0) 4 can be eliminated in favor of a term b\T^ v V \i4>S/ 'v4>, using the transfor- 



mation (3.23) parameterized by o~$. 



Alternatively, the term ciG^V »</>V v <j) can be eliminated in favor of a term biT^V «(/>V u (, 



using the transformation (3.23) parameterized by as- Our result in this representation 



agrees with that of Weinberg [18] when all the matter terms are dropped. The c\ represen- 
tation has the advantage over the b\ representation that the corrections are confined to the 
scalar sector and do not involve matter. The term ciG^V \i4>V v<f> has interesting effects: it 
can give rise to a self-tuning cosmology as well as potentially support a Vainshtein screening 
mechanism 1311. 



8 More precisely the criterion is that the zeroth order term in the expansion in a' in powers of e is nonzero. A 
nonzero a' that is proportional to e would be insufficient to allow this transformation. 
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As discussed in Appendix [Bj it is possible to eliminate all the stress-energy terms from the 
action by applying field redefinitions. This yields a form of the theory in which the weak 
equivalence principle is manifest. However, the resulting action contains higher derivative 
terms, unlike all the representations discussed so far in this subsection. As discussed in 
the Introduction and in Appendix [Cj to define the theory when higher derivative terms are 
present we use the reduction of order technique applied to the equations of motion. 

Finally, the result can be cast in the Jordan conformal frame by doing a conformal trans- 
formation, followed by some field redefinitions to simplify the answer. The result is similar 



in form to the Einstein frame action (4.5): 

"1 



(Fx, 



-g 



2 mle- a R- 2 (V4>) 2 -U(4>) 



+ e / d 4 x^g\ ~ ai {V^) A + b 2 T(V4>) 2 + ciG^V^V^ 



+ d 3 [r 2 - AR^R^ + R^pR^o) + d^^C^C. 
+ eiTT™ + e 2 T 2 



Xpa/3 



(4.6) 



Here g^ v = e a g^ u and the field is a function of <j), where the function is chosen to 
give canonical normalization for (j> in the Jordan frame action (4.6). All of the functions 
U, &i, etc. in this action differ from the corresponding functions in the Einstein frame 



representation (4.5), but can in principle be computed in terms of them. The Jordan frame 



result (4.6) can also be cast in a number of different forms using linearized field redefinitions, 



just as for the Einstein frame result (4.5). Note that the stress energy tensor we use is the 



same in both frames, and is defined in Appendix [Aj The result (4.6) matches that found 
by Park et al. [21] (up to some minor adjustments, see Appendix |P]) . 

We note that the Chern-Simons term (d^) gives rise to third-order derivatives in the equations 
of motion [see Eqs. ( E.4 ) and ( E.5 ) below] . However, with the choice of foliation^] given by surfaces 
of constant </>, there are no third-order time derivatives, and so the Chern-Simons term is not 
a higher derivative term according to our definition (see the discussion in Sec. 1.2 above), and 



is not subject to the Ostrogradski instability. For further discussion of the Chern-Simons term 
in gravitational theories, see, e.g., Ref. [32]. As a parity- violating term, this term modifies the 
propagation speed of different polarizations of gravitons. 

In the above derivation, we eliminated higher derivative terms using field redefinitions. As 
discussed in Appendix [Cj an alternative but equivalent procedure is to derive a form of the action 
which explicitly exhibits the extra degrees of freedom associated with the higher derivative terms, 
and then integrate out those degrees of freedom at tree level. This is shown explicitly for higher 
derivatives of the scalar field in Appendix [Cj and can also be shown explicitly for the terms d\ and 
d 2 involving higher derivatives of the metric. A third, equivalent method is to perform a reduction 



9 This choice requires the assumption that \7(f> is timelike everywhere, which will be true in cosmological appli- 
cations when perturbations are sufficiently small. 
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of order procedure at the level of the equations of motion, as discussed in the Introduction and 
in Appendix [Cj 

The above derivation confirms the general argument made in the Introduction that it should 
not be necessary to include higher derivative terms in the action. This is because the new terms 
that are generated when one eliminates the higher derivative terms should already be included 
in the derivative expansion. In the above derivation, if we eliminate from the start the higher 
derivative terms (03, 63, 64, 65, be, C2, C3, d\, cfe), then we must also forbid all transformations that 



generate these terms, which includes all the transformations we have considered except Eqs. (3.3), 



(3.4), (3.11) and (3.23). The above derivation gets modified by dropping steps 2, 3, and 4, the 



portion of step 5 that sets to zero, and the portion of step 6 that sets 63 to zero. The final 



result (4.5) is unchanged. 



In a similar vein, the correct result can also be obtained by omitting from the initial action all 
terms involving the stress energy tensor, that is, the terms parameterized by b\, . . . , 67 and e\, 
e-i. If one follows all the steps of the derivation in Table [2j the same final result is obtained, and 
all the final coefficients are nonzero in general. This occurs because all the terms involving the 
stress energy tensor have alternative representations not involving it (although they do involve 
higher derivatives). Thus, from this point of view, it is not necessary to include in the action 
stress-energy terms. 

However, it is not possible to do without both the higher derivative terms and the stress-energy 



terms. Suppose we throw out at the start all the higher derivative terms in both the action (2.1 ) 



and Jordan frame metric (2.4), and in addition omit all the stress-energy terms in the action. 



This would yield a version of the action (2.1) involving only the terms ai, 02, fa, fa, ci, d% and 



g?4. Using the transformation (3.2) the terms fa and fa can be exchanged for b\ and b 2 , and the 



terms 02 and b\ can then be eliminated using the transformations parameterized by 03 and o~g. 



This yields our final action (4.5) but without the terms e± and e 2 - Therefore, for a fully general 



theory, one can choose to eliminate higher derivative terms, or stress-energy terms, but not both. 



4.3 Extension to N scalar fields: Qui-N-tessence 

The preceding analysis can be generalized straightforwardly to the case of N scalar fields, which 
we call "qui-N-tessence" , an analog of multifield inflation 
gets replaced by a general nonlinear sigma model: 

,,2 

Sn = I d 4 x 



16] . The zeroth order action (2.2) 



-9 



■R 



1 



<Mb(</> C )V^V 



U(cP C ) 



(4.7) 



where (f> A = (cj) 1 , . . . , <p N ) are the N scalar fields and qAB is a metric on the target space. In the 
remainder of the action, functions of 4> are replaced by functions of 4> A . The first three terms of 



the second line of Eq. (4.5) are replaced by 



axABCDV^V^Vx^V^g^g^ + a 2 ABC V M V A </> A V ff <£ B V \ct> C 'sT 
+ c 1AB G^V,4> A V u 4> B - 



Thus the coefficients a%, a 2 and c\ become tensors on the target space of the indicated orders. 
Note that we must use the representation involving the coefficients a\ abc d > 02 ABC an d not 
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biAB,b2AB (we assume a a ^ 0) since the latter are less general; the equivalence between the 



different representations discussed in Sec. |4.2| does not generalize to the N field case. When 
N > 4 one could also add a term 

aiABCDV^V^Vx^V^e^, (4.9) 

where a^ABCD is an arbitrary 4-form on the target space. 

5 Order of Magnitude Estimates and Domain of Validity 

In the previous sections, we started from the standard quintessence model with a matter cou- 



pling (2.2), and added arbitrary corrections involving the scalar field and metric in a derivative 
expansion up to four derivatives. We then exploited the field-redefinition freedom to eliminate all 
terms that give rise to additional degrees of freedom ( "higher derivative terms" ) , and to reduce 



the set of operators in the action to the canonical and unique set given in our final action (4.5). 

We now turn to estimating the scaling of the coefficients in the final action using effective 
field theory. We will then use these estimates to determine the domain of validity of the theory. 

5.1 Derivation of Scaling of Coefficients 

We start by recalling the scenario of pseudo Nambu-Goldstone Bosons (PNGBs) [271 128] discussed 



in the Introduction that may give rise to the zeroth order action (2.2). Suppose that at some 
high energy scale M* we spontaneously break a set of continuous global symmetries and thereby 
generate N massless Goldstone bosons <j> = ((f) 1 , . . . , (f> N ). The theory then has N residual 
continuous symmetries. If we now suppose that these residual symmetries are explicitly broken 
at some much lower energy scale A, then a potential is generated that scales as 

A 4 F(0 A /M*), (5.1) 

for some function V which is of order unity. In particular the mass of the PNGB fields scale as 
A 2 /M* and can be very light. For example, in axion models M* ~ 10 12 GeV and A ~ Aqcd ~ 100 
MeV, giving an axion mass of order 10 -5 eV. 

The leading order action for the PNGB fields coupled to gravity at low energies will be that 
of a nonlinear sigma model, 

~1 2 1 

-m n R ( 

2 p 2 



S= (fx 



-m 2 p R - -q A B^ C '/MOV^V^V* - A 4 V((/> C /M* 



(5.2) 



where qAB is a metric on the target space which admits N Killing vectors (the residual sym- 
metries). In the special case where qAB is flat, these residual symmetries are shift symmetries 
<fi A — > <p Aj r constant. We now assume that these fields drive the cosmic acceleration, and in ad- 
dition we assume that the kinetic and potential terms are of the same order, that is, we assume 



that slow roll parameters are only modestly small. It then follows from the action \b.2\ that the 
scales of spontaneous and explicit symmetry breaking M* and A must be of ordei 



10 



M* ~ m p , A ~ a/ H^nip, (5.3) 



10 The need to use the Hubble scale today in the symmetry breaking scale A is associated with the coincidence 
problem. 
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where Ho is the Hubble parameter, so that the quintessence fields have mass ~ Hq and energy 
density ~ ui^Hq. Defining the dimensionless fields ip A = 4> A /m p allows us to rewrite the action 
as 



5 



a 



m z p R - ^-m z p q AB ^ c )V vV A V^g^ - rn^V^) 



B av 



(5.4) 



Consider now the stability of the theory (5.4) under loop corrections. The story is exactly 
the same here as in inflationary models 



[27] (aside from couplings to matter, see below). 
Computing loop corrections starting from the action ( |5.4[ ) does not lead to large corrections 
5m ^> Hq to the mass of the quintessence fields, because in the limit where the explicit symmetry 
breaking scale A = ^m p Ho goes to zero, the theory possesses exact symmetries which must 
be respected by the loop corrections. Hence the loop corrections to the potential must scale 
proportionally to H^m 2 , as for the original potential. Thus the smallness of the mass of the 
quintessence field is natural in the sense of 't Hooft. However, this is not the entire story, since 



the form (5.4) of the low energy theory imposes non trivial constraints on the physics at high 



energies, which must respect the residual symmetries. Indeed in general there is no guarantee 



that there exists a consistent high energy theory with the low energy limit (5.4). This question 



is beyond the scope of this paper: we shall simply assume that a consistent UV theory can be 
found. See Ref. [35] for an example of an attempt to address this issue. 

So far in the discussion we have neglected coupling to matter. If we assume the validity of 
the weak equivalence principle, the general leading order coupling of <\f to matter will be of the 



form of a scalar-tensor theory, given by adding to the action (5.4) the term 



= a(07M,) 



(5.5) 



for some function a. 

We now suppose that one or more of the PNGB fields has a mass ~ M which is parametri- 
cally larger than Hq, and we integrate out these heavier fields, following the similar analysis of 
inflationary models by Burgess, Lee and Trott [19]. Integrating out the heavier fields gives rise 
to modifications to the target space metric and potential for the remaining light fields [that do 



not change the scalings shown in Eq. (5.4)], and also a set of correction terms to the leading 



order action (5.4). The leading, tree- level correction terms can be obtained simply by solving the 



classical equations of motion for the heavy fields in an adiabatic approximation and substituting 
back into the action. One finds that the induced correction terms have the fornFH 



M d 



(5.6) 



where the sum is over operators O n involving d derivatives acting on k powers of the dimensionless 
fields <p and/or g nu , and the coefficients c n are of order unity (see Appendix |F| for details). In 
other words, each additional derivative is suppressed by a power of the mass M of the fields that 
have been integrated out (which we can think of as a cutoff scale), and the overall prefactor is 
such that the normal kinetic terms would be reproduced for the case k = d = 2. 



11 These are the terms involving just the scalar field and metric. One also finds correction terms involving the 
matter stress energy tensor as long as a' 7^ 0, of the form indicated in Table |3] 
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Coefficient 


Term in Action 


Scaling 


a i 


(V0) 4 


~ l/(m 2 M 2 ) 




USiVcb) 2 


~ l/(m„M 2 ) 
/ \ y ) 


03 t 


(Dcj)) 2 

\ T J 




h 




~ l/(m 2 M 2 ) 

1 \ y j 




T(Vcj)) 2 


~ l/(m 2 M 2 ) 

1 \ y j 


b-i t 






64 t 






h t 






b 6 t 


i?T 




b 7 


T 




Cl 




~ 1/M 2 


C2 t 


R{Vct>) 2 




C3 t 


RU(j) 




di t 


R 2 




d 2 t 


R^ v R^iu 






Gauss-Bonnet 


~ m 2 } /M 2 


C?4 


Chern-Simons 


~ m 2 } /M 2 


ei 




~ l/(m 2 M 2 ) 


e2 


T 2 


- l/(m 2 M 2 ) 



Table 3: This table gives the scalings of the various coefficients. The first column lists the coefficients, and 



the second column lists the corresponding terms in the action (2.3 1. Daggers in the first column indicate 



higher derivative terms. The third column gives our estimate of the scale of the coefficients, under the 



assumptions discussed in the text, for those coefficients that are nonzero in our final action (4.5 1, or in 



versions of that action obtained using the field redefinitions (3.11 1 or (3.23). The quantity M is the mass 



of the lightest field that is integrated out to produce our final action. In all cases, these scales for the 
coefficients correspond to fractional corrections to the leading order dynamics of order ~ Hq /M 2 . 
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Note that the rule (5.6) for how the coefficients of additional corrections to the action depend 



on the cutoff scale M differs from the usual rule of effective field theory, where an operator 

M~ D . 



The rule (5.6) instead gives a coefficient 



of dimension D + 4 has a coefficient 

M~( d ~ 2 ^mp^ k 2 \ where d is the number of derivatives in the operator and k is the number of 
powers of (canonically normalized) fields, related to D by D = d + k — 4. The difference between 
the two rules arises from the fact that we are making nontrivial assumptions about the physics 
above the scale M, specifically that it is described by an action of the PNGB form (5.4) 12 If we 



were to allow arbitrary physics at energies above the scale M, then the coefficients would scale 
according to the standard rule. 



We now specialize to the case of a single light field. The correction terms (5.6) have the form 



of a double power series, in number of derivatives and in powers of the fields. If we fix the number 
of derivatives and associated index structure, we can sum over all operators that differ only by 
powers of cp to obtain operators with prefactors that are functions of tp, f(<p) = Ck<p k with 
coefficients of order unity. We now write out all the resulting terms to leading order in 1/M 2 , 



imposing general covariance. The result is the theory (4.5) discussed in the last sectior 13 , but 



with additional information about the coefficients oi, b\ etc. Specifically we find that 



1 



ail 



m? p M 2 



aiO/rap), 



(5.7) 



where the function di is of order unity, i.e., the coefficients in its Taylor expansion are independent 
of m p and M. The corresponding prefactors or overall scaling for the other coefficients are listed 
in Table [HI 

Finally, we note that, as is well known, Solar System tests of general relativity strongly 
constrain the coupling of <fi to the matter sectoi 1 

A = m p \a'(4>o)\, where <f>Q is the present day cosmological background value of 
System constraint ra^lA < 10 -2 [30]. In addition the coupling of the scalar to the visible sector 



If we define the dimensionless parameter 

then the Solar 



will generically give rise to large corrections to the quintessence potential via loop corrections 
|41| I32"t |4"51 14"4"1 14"5] . For a fermion of mass m / , the correction 5m to the mass of the quintessence 
field will be of order 



5m 



J m2 f ) 
yH m p J 



(5- 



More general interactions which are not of the form \5.4\ can modify the scaling rule (5.6 1, even if they respect 
the residual (shift) symmetries. For example consider a scalar field ip of mass m which couples to <j> via a term 
ip(V '(f>) 2 /m„ for some mass scale m„. Integrating out this field gives a correction to the (f> action ~ (V(/>) 4 /(m 2 m 2 ) 
(see Appendix [fJ . To keep such terms from invalidating the scaling rule we need to assume that mm* > Mm p , 
i.e. that any such fields are either sufficiently massive or sufficiently weakly coupled to the PNGB fields. 

13 The parity-violating Chern-Simons term is not generated in this way, since the fields we are integrating out do 
not violate parity. To obtain the Chern-Simons term with the scaling indicated in Table[3]would require integrating 
out some parity violating fields at the scale M which approximately respect the residual (shift) symmetries. 

14 Strictly speaking, Solar System tests lie outside the domain of validity of our effective field theory unless 
M" 1 < 1A.U., which is very small compared to H ; see Sec. 

are important in the Solar System, such 



above. 

15 This constraint can be evaded in models where nonlinear effects in 
as Chameleon [36] and Galileon models 37, [3811391. 
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If A ~ 1 and mf ^> ^m p Ho ~ eV, then 5m 3> Hq, which is inconsistent if the quintessence 
field is to drive cosmic acceleration. This is a well known naturalness problem for matter couplings 
in quintessence models, and it motivates setting^] a = 0. 



5.2 Domain of Validity of the Effective Field Theory 



We now estimate the domain of validity for the theory (4.5) with the scalings given by Table 
[3j by requiring that the terms with higher derivatives be small compared to terms with fewer 
derivatives. If E is the energy involved in a given process, or equivalently E~ l is the corresponding 
time-scale or length-scale, then successive terms in the derivative expansion are suppressed by 
the ratio E/M, which yields the standard condition 

E<^M (5.9) 

for the domain of validity. As discussed in the Introduction, M must be somewhat larger than 
Hq in order to describe the background cosmology and observable perturbation modes. However 
if M is significantly larger than Hq then the corrections due to the higher order terms in Eq. 



(4.5) become negligible, and the theory reduces to a standard quintessence model with some 
matter coupling. Therefore, the interesting regime is when M is perhaps just one or two orders 
of magnitude larger than Hq, as indicated in Fig. [TJ In particular, when the scale M is in this 
interesting regime, the theory is unable to describe gravitational effects in the Solar System and 
binary pulsars, which is a shortcoming of the effective field theory approach used here. 



Consider now the background cosmological solution. The theory (2.1 ) to zeroth order in e (or 
equivalently 1/M 2 ) has the equations of motion 

m 2 p G al3 = V a (jN P (t> - \{V4>?g a p ~ U(4>)g a p + e 2a ^T a p, (5.10a) 
U(j) = U'{</>) - ^a'e 2a T. (5.10b) 

For each of these two equations we assume that all of the terms are of the same order. For the 
matter terms this is this is a reasonable approximation, since f^A ~ 0.7 and ^matter ~ 0.3. If 
the scalar potential term dominates over the kinetic term, then the following estimates need to 
be modified by including factors of slow roll parameters; we ignore these factors here since we 
expect them to be only modestly small. Similarly, our estimates assume that m p a' is of order 
unity; some changes would be required if this quantity were very small. From these assumptions, 
and ignoring 0(1) functions of the scalar field, we have 

m 2 p R ~ (V</>) 2 ~ U ~ T ~ m p U§ ~ m p U'(4>) ~ Hlm 2 p . (5.11) 



Inserting these estimates into the action (4.5) and using the scalings given in Table |3| we find 
that for each of the correction terms in the action, the fractional corrections to the leading order 
cosmological dynamics scale as Hq/M 2 . The corrections therefore are of order unity at M ~ Ho, 
as we would expect, since at this scale the heavy fields which we have integrated out have the 



16 More precisely the condition is a' = 0, i.e., a — constant, but the constant can be absorbed by a rescaling of 
all the dimensionful parameters in the matter action. 
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same mass scale as the light fields, and would be expected to give rise to 0(1) corrections to the 
dynamics. This gives a useful consistency check of the calculations underlying Table [3] discussed 
in the previous subsection. 



In addition to the standard constraint (5.9), there are other constraints on the domain of 



validity which we now discuss. We focus attention on cosmological perturbations, for which 
4>(t, x) = (j>o (t) + S(p(t, x), and consider the conditions under which the dynamics of the perturba- 
tion 5(p can be described by the effective theory. Consider localized wavepacket modes 5(p, where 
the size of the wavepacket is of the same order as the wavelength, both ~ E . For such modes 
we can characterize perturbations in terms of two parameters, the energy E and the number 
of quanta or mode occupation number N . The total energy of the wavepacket will be of order 



NE ~ J d 3 x(V5(f)) 2 ~ E 3 (E5(f>) 2 which gives the estimate 



NE. 

The fractional density perturbation due to the wavepacket is of order 

{V5(f>) 2 NE A 



Sp 
P 



H o m P 



2 ' 



(5.12) 



(5.13) 



We now demand that the term ai(V<5^) 4 in the actior 17 be small compared to the leading order 
term (VScp) 2 . Using the scaling a\ ~ l/{m 2 M 2 ) from Table [3] and combining with the estimate 



(5.13) of the fractional density perturbation then gives the constraint 18 

Sp 



P 



(5.14) 



Thus, the theory can describe perturbations in the nonlinear regime, but the perturbations can 
only be modestly nonlinear if M is fairly close to Hq. In terms of the parameters E and N the 



constraint (5.14) is 



NE 4 < M 2 ml. 



(5.15) 



This gives a nontrivial constraint on the domain of validity of the theory in the regime E < M. 
The two dimensional parameter space (E,N) is illustrated in Fig.[2j which shows the constraints 



(5.9) and (5.14), the curves Sp/p ~ 1 and bp ~ M 2 /Hq, as well as the curve where 



~ m r . 



Another potential constraint on the domain of the validity of the theory (4.5) with the scalings 
given by Table [3] is that the theory should be weakly coupled, i.e. the effects of loop corrections 
should be small. Using the power counting methods of Ref. [19] one can show that this is indeed 



17 Here we envisage computing an action for the perturbations by expanding the action (4.5 I around the back- 
ground cosmological solution, as in Ref. [19!. 

18 In the previous subsection we showed that ai(4>) — a\ ((f>/m p )/ (M 2 m P ), where fii is function for which all the 
Taylor expansion coefficients are of order unity. It follows that a\ 



1 for 



m p . However the estimate (5.141 



requires the stronger assumption a\ < 1 for (j> 3> m p which need not be valid. If we instead assume that a\ ~ 
(<j}/m p ) a for <f> » m p then the constraint |5.15| gets replaced by N(E/M) 1 < m 2 p /M 2 , where 7 = 2 (4 + a)/ (2 + a). 



This modifies the boundary of the domain of validity of the effective field theory shown in Fig. [2] by changing the 
slope of the tilted portion of the boundary. In the limit a — > 00 this portion of the boundary approaches the green 
curve Sip ~ m p . 
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true within the domain Hq < E <C M of interest. Strong coupling can arise due to tri-linear 
couplings, as discussed in Sec. 2.2 of Ref. [19], but this only occurs for energies far below the 
Hubble scale Hq, and so is not relevant to cosmological applications of the theory. 

We note that there are several well known theories of cosmic acceleration that are not en- 
compassed by our effective field theory. The form of our expansion requires that the dominant 
contribution to cosmic acceleration be the leading order scalar terms and not the higher order 
terms, and so theories in which other mechanisms provide the acceleration cannot be described 
in our formalism. One example is provided by K-essence models in which terms in the action 
like (V</>) 4 , (V0) 6 . . .are all equally important. In particular this is true for ghost condensate 
models [6]. Also there are many cosmic acceleration models that exploit the Vainshtein effect 
[HI 27HH] t° evade Solar System constraints on light fields with gravitational-strength couplings. 
The Vainshtein effect relies on nonlinear derivative terms in the scalar field action. Although our 
class of theories includes models that demonstrate the Vainshtein mechanism, the mechanism 
only operates outside the domain of validity of our approach, as we require the nonlinear deriva- 
tive terms to be small. The chameleon mechanism [S6\ 149] . on the other hand, does not require 
nonlinearities in the derivatives of the scalar field, and thus may be analyzed in our formalism, 
although the regime in which a screening mechanism would be required to evade fifth force ex- 
periments and solar system constraints will be in the regime of validity of our analysis only for 
large enough values of the cutoff M. 

6 Discussion and Conclusions 

In this paper, we have investigated effective field theory models of cosmic acceleration involving 
a metric and a single scalar field. The set of theories we considered consists of a standard 
quintessence model with matter coupling, together with a general covariant derivative expansion, 
truncated at four derivatives. We showed that this class of theories can be obtained from a PNGB 
scenario, where one of the PNGB fields is lighter than all the others, and the heavier fields are 
integrated out. We showed that in constructing this class of theories, including higher derivative 
terms in the action, as suggested by Weinberg [18], does not give any increased generality. We 
also showed that complete generality requires one to include terms in the action that depend on 
the stress-energy tensor of the matter fields. 

We now turn to a discussion of some of the advantages and shortcomings of the approach 
adopted here to describe models of dark energy. Some of the shortcomings are: 

• By construction, our approach excludes theories where nonlinear kinetic terms in the action 
give an order unity contribution to the dynamics, such as K-essence, ghost condensates etc., 
since such theories do not arise from the PNGB construction used in this paper. On the 
other hand, such theories are less natural than the class of theories considered here, from 
the point of view of loop corrections: they require very nontrivial physics at the scale ~ Hq, 
instead of at the scale ~ ^HQm p required in the PNGB approach. The most general class 
of theories of this kind is that of Horndeski |50] , which contains four free functions of <j) and 
(V</>) 2 |13] , and which is the most general class of theories of a metric and a scalar field for 
which the equations of motion are second order. As discussed in the Introduction, these 
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theories are included in the alternative, background-dependent approach to effective field 
theories of quintessence of Creminelli et al. |20j. 

• Our class of theories will be observationally distinguishable from vanilla quintessence theo- 
ries only if the cutoff M is near the Hubble scale Hq. In this regime, our framework cannot 
be used to analyze Solar System tests of general relativity, since they are outside the domain 
of validity of the effective field theory. Also, when the background cosmology is evolved 
backwards in time it passes outside the domain of validity at fairly low redshifts. (This is 
not a serious disadvantage since dark energy dominates only at low redshifts.) 

• We have restricted attention to theories with a metric and a single scalar field, with the 
only symmetry being general covariance. Thus, our analysis does not include models with 
several scalar fields, vector fields etc. In addition, our analysis excludes an interesting class 
of models that one obtains by imposing that the action be invariant under <j) — > /(</>), where 
/ is any monotonic function, as such a symmetry cannot be realized with our derivative 
expansion. This class of models includes Horava-Lifshitz gravity and has the same number 
of physical degrees of freedom as general relativity |13| 51J. It would be interesting to 
explore the most general dark energy models of this kind. 

Some of the advantages of the approach used here are: 

• Our class of theories is generic within the PNGB construction, which itself is a well moti- 
vated way to obtain the ultralight fields needed for cosmic acceleration. The theories are 
fairly simple and it should be straightforward to confront them with observational data. 

• Our class of theories allow for a unified treatment of the cosmological background and 
perturbations, unlike the background-dependent approach of Ref. [20] , 

Finally, we list some possible directions in which the approach used here could be extended: 

• It would be interesting to compute the relation between the nine free functions used in our 
theories to the free functions of the post-Friedmannian approach to parameterizing dark 
energy models [15] . 

• It would be interesting to explore the phenomenology of the various higher order terms 
in our action, for the cosmological background evolution and perturbations. Many of the 
terms have already been explored in detail, see for example Refs. [30] 131] . 

• Either by using the post-Friedmannian approach, or more directly, it would be useful to 
compute the current observational constraints on the free functions in the action. 

• An interesting open question is the extent to which our final action is generic. That is, is 
there a class of theories more general than nonlinear sigma model PNGB theories for which 
our action is obtained by integrating out some of the fields? 
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A Notation and Conventions 



We use natural units with c = h = 1, and define the the reduced Planck mass via m 2 = 1/8ttG. 
The Einstein and Jordan frame metrics are and g^ u respectively, and the corresponding 
derivative operators are V M and V M . We use the usual abbreviations (V^) 2 = g^V M </>V ' v i/j and 
□0 = V^V^. Primes denote derivatives with respect to the scalar field eft, as in U'(4>). We use 
the (— , +, +, +) metric signature and the sign conventions (+, +, +) in the notation of Ref. [52] . 
Finally we take e* 1 p to be the antisymmetric tensor with g 0123 = \ j yj— g. 

We define the (Jordan- frame) stress-energy tensor T^ v in the usual way in terms of the 
Jordan- frame metric g^ that appears in the matter action S* m : 

S m [g„u + Sfa, ^n] - S m [g^, M = \J d 4 x^T^g» x 5g Xu + 0(Sg 2 ). (A.l) 

We then define T = T^, and define the quantities T^ v and T^ u by raising and lowering indices 



with the Einstein-frame metric g^u, which is related to g^ u via Eq. (2.4). To zeroth order in e 
this stress energy tensor obeys the conservation law 

e- 2a V x (e 2a T Xa ) = X -dTV a <\> + O(e). (A.2) 

B The Weak Equivalence Principle 

In this Appendix, we show that including terms in the action that depend explicitly on the 



matter stress energy tensor, as in Eq. (2.1) above, generically gives rise to violations of the weak 



equivalence principle. However, we also show that our specific model (2.1) does not, to linear 
order in e. Since the parameter e essentially counts the number of derivatives in our derivative 
expansion, it follows the weak equivalence principle is satisfied for our derivative expansion up 
to four derivatives. 

B.l Generic Violations of Weak Equivalence Principle when Stress-Energy 
Terms are Present in Action 

Consider first an action principle of the general form 

S[g a pA,^m\ = S % [g a p,(f>] + S m [gap, ip m ]. (B.l) 

Here the first term is a gravitational action, depending only on the metric g a p and the scalar 
field (f>, and the second term is the matter action, in which all the matter fields ip m couple only 
to the Jordan metric g a p (some function of g a p and <j>), and not to g a p and cf> individually. By 
definition, any theory of this form obeys the weak equivalence principle. What this means is 
as follows. We define weakly self-gravitating bodies to be bodies for which we can neglect the 



perturbations they cause to g a p and <j). From the form of the action (B.l), it follows that all 
weakly self-gravitating bodies will fall on geodesies of the metric g a /3, and hence will all fall on 
the same geodesies. 



The action principle (2.1) we use in this paper is not of the general form (B.l), because of 



the explicit appearance of terms involving the stress energy tensor in the gravitational action. 
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Therefore one expects violation of the weak equivalence principle to arise. We now verify explicitly 
that this does occur in a specific example. We choose the following special case of the action 



(2.1), where the only perturbative term included is the term proportional to the trace of the 



stress energy tensor: 
S~- 



-9 



1 



-miR 



1 



(V^) 2 - U(<t>) + ef(4>)T 



2 p 2 

We choose the matter field ip m to be a scalar field ij) with action 

"1 



+ S m [g a p,i/) 



d*Xy 



-9 



(B.2) 



(B.3) 



9: 



and we specialize the relation (2.4) between the two metrics to be the conformal transformation 
- e a Wg a g. This gives T = — e _a (V?/>) 2 — 4V and the total action is therefore 



/ 



d A XyJ—g 



1 



-miR 



1 



1 



(V<PY - U{<t>) --{e a + 2ee- a f){V^Y ~ (e 2a + 4e/)Fty>) 



(B.4) 



2" p - 2 V ' T ' ~ VT/ 2' 

The kinetic term for tp can be written as J drxy/^g(VipY where g a p = (e a + 2ee~ a f)g a /3, and 
the potential term can be written as J d^x^—gVfy), where g a/ 3 = \J e 2a + ^efg a /3. Therefore, 
objects whose stress energy is composed of different combinations of the kinetic term and the 
potential term will fall on different combinations of the metrics g a/ s and g a R, violating the weak 
equivalence principle. 

B.2 Validity of Weak Equivalence Principle to Linear Order 

In the above analysis, we note that the metrics g a n and g a p coincide to linear order in e, so there 
is no violation to this order. We now show that, similarly, none of the stress-energy-dependent 



terms included in Eq. (2.1) violate the weak equivalence principle, to linear order in e. 



The key idea of the proof is to use the transformation laws derived in Sec. [3] above to rewrite 



the theory in the general form (B.l), which we know satisfies the weak equivalence principle. 



All of the terms in the action given by Eqs. (2.1) - (2.4) are of this form, except for the terms 
parameterized by the coefficients &i, . . . , 67, e\ and However, as we now show, we can use 
transformations to eliminate these terms in favor of the remaining terms which manifestly satisfy 
the principle. 



Consider first the terms in the action (2.3) which depend linearly on the stress-energy tensor. 



We can eliminate the terms parameterized by 61,..., b$ using the transformation (3.3) with 
Pi = —2e~ 2a bi for 1 < i < 6. This generates contributions to the the terms parameterized by 
fix, . . . , /3q in the definition (2.4) of the Jordan metric. Similarly, by using the transformation 



(3.4) with a = -2e~ Ia b 7 , we can eliminate the term parameterized by 67 in favor of an 0(e) 



correction to the function a in Eq. (2.4). 



We now turn to the terms in the action (2.3) which depend quadratically on the stress-energy 



tensor, namely the terms parameterized by e\ and eg. For e\ we use the transformation (3.29) 
-e~ 2a ei, and for e2 we use the transformation (3.27) with uio = ■ 



with o"ix = 

transformations generates new contributions to the linear stress-energy terms parameterized by 



201 e 2 . These 
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bi, 62 > ^5) be and 67 (see Table [T]), but we have already shown that all of those terms satisfy the 
weak equivalence principle. 



To summarize, we have shown that our model (2.1) satisfies the weak equivalence principle 
despite the explicit appearance of stress energy terms in the action. Of course, there can be 
violations of the strong equivalence principle in models of this kind, which can even be of order 
unity [53j. In addition, the weak equivalence principle will generically be violated by quantum 
loop corrections, although this is a small effect 



B.3 Potential Ambiguity in Definition of Weak Equivalence Principle 

We next discuss a potential ambiguity that arises in the definition of the weak equivalence princi- 
ple. In the definition one restricts attention to bodies whose gravitational fields, as measured by 
the perturbations they produce to the metric g^ v and scalar field <j>, can be neglected. However, 



consider for example the field redefinition (3.27), where the metric transforms according to 



9al3 = 9a/3 + ^ w Tg a p. (B.5) 
It is possible for the perturbation 5g a R generated by the body to be negligible, but the perturba- 



tion 5g Q p to be non-negligible, because of the appearance of the stress-energy term in Eq. (B.5). 
If this occurs then the weak equivalence principle could be valid for one choice of variables, but 
not valid for the other choice. 

To assess this ambiguity, we now make some order of magnitude estimates. Consider a body 
of mass ~ Mb and size ~ R. Then in general relativity the size of the metric perturbation due 
to the body is of order 5g a p ~ Mb/irriiR). Suppose now that g\q ~ l/{m 2 M 2 ), as indicated by 
Eq. (3.28) and Table [3j Then the contribution to the metric perturbation 5g a p from the second 
term in Eq. (B.5) will be of order M&/ (R 3 m 2 M 2 ) , which will be much larger than 5g a p whenever 



R <C M . Therefore the ambiguity could in principle arise. 

However, in the models considered in this paper the ambiguity does not occur. This is because 



the condition R <C M 1 is excluded by the condition (5.9) for the validity of the effective field 
theory. 

C Equivalence Between Field Redefinitions, Integrating Out New 
Degrees of Freedom, and Reduction of Order 



The action (2.1) we start with in the body of the paper contains several higher derivative terms, 
that is, terms which gives contributions to the equations of motion which involve third-order 
and fourth-order time derivatives of the fields. As discussed in the Introduction, the theory with 
these higher derivative terms contains additional degrees of freedom compared to our zeroth order 



action (2.2), which contains a single graviton and scalar. In this paper our goal is to describe a 



general class of theories containing just one tensor and one scalar degree of freedom, so we wish 



to exclude these additional degrees of freedou 19 



19 Higher derivative terms are also generically associated with instabilities [22], although this can be evaded in 
special cases, for example R 2 terms. 
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Therefore, as discussed in the Introduction, we define the theory we wish to consider, associ- 
ated with our action (2.1), to be that obtained from the following series of steps: 



1. Vary the action to obtain the equations of motion, which will contain third-order and 
fourth-order derivative terms which are proportional to e. 

2. Perform a reduction of order procedure on the equations of motion |23} |2~H I25j . That is, 
substitute the zeroth-order in e equations of motion into the higher derivative terms in order 
to obtain equations that contain only second-order and lower order time derivatives, which 
are equivalent to the original equations up to correction terms of 0(e 2 ) which we neglect. 

3. Optionally, one can then derive the action principle that gives the reduced-order equations 
of motion. 

In this Appendix, we show that this procedure is equivalent to the computational procedure we 
use in the body of the paper, in which we apply perturbative field redefinitions directly to the 
action in order to obtain an action with no higher derivative terms. We also show that it is 
equivalent to integrating out at tree level the extra degrees of freedom that are associated with 
the higher derivative terms. 

We note that the analyses of general quintessence models by Weinberg [18] and Park et al. 
|21j used a different method of eliminating higher derivative terms. They performed a reduction 
of order procedure directly at the level of the action, that is, they substituted the zeroth-order 
equations of motion directly into the higher derivative terms in the action, to obtain an action 
with no higher derivative terms. We will show that this method is not in general correct; it does 
not agree with the theory obtained by applying the reduction of order method to the equations 
of motior 20 However, it differs from the correct result only by field redefinitions (that do not 
involve higher derivatives), and so for the purpose of attempting to classify general theories of 
quintessence, Weinberg's method is adequate. 

C.l Reduction of Order Method 

We start by considering the case of just a scalar field; a more general argument valid for scalar 
and tensor fields will be given below. Consider a general action of the form 

S = j d 4 x^-g j-^) 2 " U(4>) + eF[0,(V0) 2 , □<£]}, (C- 1 ) 

where F is an arbitrary function. We introduce the notation K = (V0) 2 and L = \3(p. We first 
show that applying the reduction of order procedure to the equations of motion (steps 1-3 



above) give rise to a theory of the form (C.l ) but with F((p, K, L) replaced by another function 
F(<f>,K,L), given by 

F(4>, K, L) = F[4>, K, U' (</>)] + [L- U'(<l>)]F L [<t>, K, U' (</>)} . (C.2) 



20 The reason is that substituting the zeroth order equations of motion into the action gives an action which is 
correct off-shell to O(e ) and on-shell to O(e), but it needs to be valid off-shell to 0{e). 
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To see this, we vary the action (C.l) to obtain the equation of motion 



D<f> - U'(<f>) + eFj - 2eV a (F K V a <f>) + eUF L = 0. 
We now make the field redefinition 

V> = + eFJ0, (V</>) 2 , □</>]. 



Rewriting the equation of motion (C.3) in terms of ip yields 



Dil> - U'{iP) + eU"[$)F iL + eF )<jt - 2eV Q (F^V a V) = 0(e 2 
where the arguments of F^,F^l and F : k are now [tp, (yip) 2 , dip]. 



(C.3) 



(C.4) 



(C.5) 



We now apply the reduction of order procedure to the equation of motion given by Eqs. (C.4) 



and (C.5), that is, we substitute in the zeroth order equation of motion Otp = U'(ip). The field 



redefinition (C.4) gets replaced by the following field redefinition which does not involve higher 
derivatives: 



V> = <P + eFjXfr (V<P) 2 , U'{4>)\ + 0(e 2 



(C.6) 



The equation of motion (C.5) is unchanged, except that the arguments of F^F^ and F ; k are 
now [ip, (S/tp) 2 , U'(ip)]. This equation of motion can be obtained from the action 



S 



d^xyj—g 



1 



(W) 2 - W) + eFty, (VV) 2 , U'm 



(C.7) 



Finally we rewrite this action in terms of <j> using the change of variable (C.6). The result is 
1 



S 



d^x^ 



-g 



-{Vcpf - U{4>) + eF[<t>, (V(/>) 2 , U'(<t>)] + ep0 - U' 



(V0) 2 ,C/' 



(C.8) 



Note that although this action contains second order derivatives, the corresponding equations of 
motion contain derivatives only up to second order, that is, the theory is no longer a "higher 



derivative" theory |30j. The final, reduced-order action (C.8) is of the form (C.2) claimed above. 



The final result (C.8) shows explicitly that the method of reducing order directly in the action 



used in Refs. p2J[2T] is not correct. Applying this procedure to the action (C.l) would yield the 



first three terms in the action (C.8), but not the fourth term. 



C.2 Method of Integrating Out the Additional Fields 



We next show that the same result (C.8) can be obtained by integrating out the new degrees of 



freedom that are associated with the higher derivative terms. Starting from the action (C.l), we 
introduce an auxiliary scalar field tp and consider the action 



s[cPM 



d 4 x, 



1 



(V<P) 2 - U{4>) + eF[cP, (V0) 2 , ^] + e(D<t> - i/j)F L [<f>, (V^) 2 , # 



(C.9) 
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The equation of motion for from this action is ip = □(/>, assuming Fj,i ^ 0, and substituting 



this back into the action (C.9) yields the action (C.l). Thus the two actions are equivalent 
classically. 

We now proceed to integrate out the field iJj, at tree level, i.e., classically. The equation of 
motion for (f> is ip = U'((f>) + O(e), and substituting this back the action (C.9) gives the same 



result (C.8) as was obtained from the reduction of order method. 



C.3 Field Redefinition Method 

We next turn to a discussion of the method we use to eliminate higher derivative terms in the 
body of the paper, using perturbative field redefinitions. That method is not generally applicable, 
but when it can be used, it is equivalent to the method of reduction of order (steps 1-3 above), 



as we now show. We start with an action of the form (C.l ), with the function F chosen to be of 
the form 



F{<j>, K, L) = g(cf>, K) + [L- U'(<j>)]h(<l>, K, L), 



(CIO) 



for some functions g and h. This is the most general form of F for which the field redefinition 
method can be used to eliminate the higher derivatives, and is sufficiently general to encompass 
the cases used in the body of the paper. First, we apply the reduction of order method. Inserting 



the formula (C.10) into Eq. (C.2) shows that the reduced-order action is characterized by the 



function F given by 



F((f>, K, L) = g(c/>, K) + [L- U\c/>)]h[(f>, K, U'{4>)\. 



(C.ll) 



However, the same result is obtained by starting with the action given by Eqs. (C.l) and (C.10) 
and performing the field redefinition 



4>^ct> + eh[<f>, (V^) 2 , U' (</>)] - eh[<P, (V^) 2 , D<f>\. 



(C.12) 



This shows the reduction of order and field redefinition methods are equivalent. 

We now give a more general and abstract argument for the equivalence, valid for any field 
content. Suppose we have a theory containing higher derivative terms in the action, proportional 
to e. Suppose that we can find a linearized field redefinition, involving higher derivatives, that has 
the effect of eliminating all higher derivative terms from the action. We can then consider this 
process in reverse: starting from a theory which is not higher derivative, by making a linearized 
field redefinition we obtain another theory which has higher derivative terms, proportional to 
e. However, the change in the action induced by the field redefinition must be proportional to 
the equations of motion. Hence, these higher derivative terms will be eliminated by applying 
Weinberg's method of substituting the zeroth order equations of motion into the 0(e) terms in 
the action. As we have discussed, Weinberg's procedure is valid up to a field redefinition of the 



type (C.6) which does not change the differential order. 
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D Comparison with Previous Work 



In this Appendix we compare our analysis and results to those of Park, Watson and Zurek |21j . 
who perform a similar computation with similar motivation, but obtain a somewhat different 
final result [Eq. (1) of their paper]. The main differences that arise are: 

• They work throughout in the Jordan frame, whereas we work in the Einstein frame. This 
is a minor difference which only affects the appearance of the computations and results, 
since it is always possible to translate from one frame to another. 

• As discussed in the Introduction and in Appendix [Cj they use Weinberg's method of elim- 
inating the higher derivative terms, consisting of substituting the zeroth order equations 
of motion into the higher derivative terms in the action, whereas we use the field redefini- 
tion method. The two methods are not equivalent for a given specific theory with specific 
coefficients, but are equivalent for the purpose of determining a general class of theories. 

• After eliminating higher derivative terms, their result is an action [Eq. (5) of their paper] 
that contains eleven functions of the scalar field, whereas our corresponding result ( |4.6| ) 
has only nine free functions. However, this is a minor difference: their function Z((p) can 
be eliminated by redefining the scalar field to attain canonical normalization, and their 



function f(4>) can be eliminated by the transformation used in step 7 in Sec. 4.1 above 



Another minor difference is that in their analysis they have in their action a Weyl squared 
term oc C a pysC a ^ ', which is unaffected by any of the transformation they make to the ac- 
tion. This Weyl squared term gives rise to higher derivative terms in the equation of motion 
that are associated with ghost-like additional degrees of freedom [55] . In our analysis the 
Weyl squared term is replaced by the Gauss-Bonnet term, which is not a higher derivative 
term, because it would be a topological term if it were not for the (^-dependent prefactor. 



Aside from the above minor differences, our result (4.5) is equivalent to the result given 
in Eq. (5) of their paper. Two major differences arise subsequently in the estimates of the 
scalings for the coefficients of the operators in the Lagrangian. 

First, Park et al. use the standard effective theory scaling rule wherein an operator of 



dimension 4 + n has a coefficient ~ Ar n , where A is the cutoff. As discussed in Sec. 5.1 
above, this corresponds to placing no restrictions on the theory that applies above the cutoff 
scale A. By contrast, our approach does place restrictions on the physics at scales above 



A, and yields the modified scaling rule (5.6). As a consequence, our cutoff A (which we 



denote by M in the body of the paper) can be taken all the way down to the Hubble scale 
Hq ~ 10~ 33 eV, whereas their cutoff must be larger than ~ ^JH§m v ~ 10 -3 eV. 

Second, Park et al. actually assume separate cutoffs for the gravitational, matter and scalar 
sectors of the theory, and estimate how each of their coefficients scale as functions of these 
three cutoffs. We do not understand completely their method of derivation of these scalings, 
but we do note that some of their scaling estimates are inconsistent with how the coefficients 
transform into one another under field redefinitions as discussed in Sec. [3] above. They then 
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proceed to drop some terms which their scalings indicate are subdominant, and arrive at a 
final action [Eq. (1) in their paper] which differs from ours, being parameterized by three 
free functions rather than nine. 



E Equations of Motion for Reduced Theory 



In this Appendix we compute the equations of motion for our final action (|4.5|), with the e\ and 
e2 terms omitted. We 
This yields the action 



e2 terms omitted. We start by using a transformation of the form (3.2) with 02 = — 2e 2a &2- 



S 



-9 



771 1 

- - (V^) 2 - U{<j>) + Gl (V(/>) 4 + cxG^V^V^ 
+d 3 (R 2 - tBTR^ + R^pR^) + d^C^Cx^] 



3 Q W^(1 + /3(V0) 2 ),^ 



(E.l) 



Here we have defined (3 = 2e~ 62; this was denoted @2 in the body of the paper. We have also 
set e = 1 for simplicity. The representation (E.l) is more convenient than (4.5) for computing 
the equations of motion since it avoids varying of the stress-energy tensor. 



Next, we vary the matter action in Eq. (E.l) using the definition (A.l) of the stress energy 
tensor T pu and the definition (2.4) of the Jordan metric g^. This yields 

5Sm = ~\J d 4 xV=ge 2a {SgT [T^ + 2T^/3(V0) 2 - pTV^V^] 

+5(p [-a'T + 2a'/3T(V(/>) 2 + (3'T(V(/)) 2 + 2/3V M TV^0 + 2pTD<f>] } . 

(E.2) 

Combining this with the variation of the gravitational action gives the equations of motion 
D<f> = U'{4>) - X -e 2a dT + 4ai [(V0) 2 D0 + 2V /1 V^V' t 0V>] + 3a[(V<i)) A + ^G^V v <f> 



+ IcxG^V^A - 4 (R 2 - 4R^R^ + - d'^C^C. 



J \paP 



+ -e 2a [2a'/3T(V0) 2 + (3'T(V(/)) 2 + 2/3V M TV M + 2/3TD0] 



(E.3) 



and 



mlG 



e 2a T, u + V M 0V„< 



-(V<t>) 2 + U{cj>) 



+ g^ciG aX V a( pV x <P - 4c 1 i? (j(/ ,V, ) </>V CT + Cl {R^{V4>? + RV^V U <P) 

- ^V CT V A (ci V ct 0VV) + ^n[ci(V</») 2 ] + 2V A V (At (ciV,)^V A 0) - V^V,[ Cl (V(/>) 2 

- □ (c 1 V M 0V ! ,</>) + 2RV i y v d z - 2g^ROd 3 + 4R„ v Dd 3 - 8R a v) V a d 3 

+ 4g^R ap V a V p d 3 + 4R p ^ a V p V a d 3 + 16C^ + 2e 2Q r M ,/3(V0) 2 - e 2a /3TV^V^. 



(E.4) 
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Here the tensor C pu comes from the Chern-Simons term, and is defined by 

= {V <I d A )e aXp ^V p R v \ + {VvVxdiYR^* (E.5) 
where *Ri lucrX = e aXpT R^ pT /2. Note that the zeroth order terms involving the stress-energy 



tensor depend implicitly on (3 through the expression for the Jordan metric given in Eq. (E.l). 

The terms involving c\ are written in the most compact manner we could find. Although it 
looks unlikely, the higher order derivatives in these terms do cancel; the full expansion of these 
terms is 

2cig pv R uX V a (tNxct> - -c x5 ^i?(V0) 2 - teyR^V <j> + Cl R pu {V 4>f + cyKV rfV v <f> 
+ 9^u [ciV ff 0V A 0V ff vV + ciV^Va^VV - ci(V0) 2 Q£ - d (□</>) 
- 2ciV (7 V M 0V fJ V^ - 2c , 1 V ff 0V (M 0V ly) V CT (/> + c / 1 V^V^(V0) 2 + c[V „<t>n<f> 
+ 2ciV M V I ^D<£ + 2 Cl V X (/)V a 4>R aflu x- (E.6) 

F Scaling of Coefficients Obtained by Integrating Out Pseudo- 
Nambu-Goldstone Fields 



In this Appendix we give some more details of the derivation discussed in Sec. 5.1 of the scaling 
of the coefficients of the operators in the Lagrangian. We divide the PNGB fields $ A into two 
groups, a set \ a with mass ~ Hq and a set ip r with mass ~ M, where M 3> Hq: 

d> A = ( X a ,V r ). (F.l) 
We assume an action for these fields of the form 

5 = | d A x^g | in: - l -q AB ^ A )V^ A V v ^ B g^ - H 2 V ( X a , } • (F.2) 



This is the same as the action (5.4) of Sec. 5.1 above, except that an extra factor has been 



inserted into the potential to make the fields have mass ~ M rather than ~ Hq, and we 
have specialized to units where m p = 1. We assume that the target space coordinates have been 
chosen so that the potential is minimized at = 0, i.e. 

V, r = (F.3) 

at V> r = 0. 

We now want to let M become large and integrate out the fields ip r at tree level. This can be 



done by using Feynman diagrams and using power counting , as in Ref. |19j . Alternatively and 



more simply, it can be done by writing out the equations of motion for the fields and invoking 



21 We note that Burgess el al. [T^ write down a scaling rule in their Eqs. (2.3) and (2.5) which is identical to 
our scaling rule (5.61 except that it is suppressed by an overall factor of M 2 jra^ for d > 2, where d is the number 



of derivatives. They say in their footnote 2 that this rule comes from integrating out a PNGB field of mass M. 
However we find that the detailed power counting calculations given in the second example in their Sec. 2.2 actually 
yield our scaling rule rather than theirs. 
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an adiabatic approximation. At zeroth order in 1/M, the theory obtained for the fields x a 1S a 



nonlinear sigma model where the potential is just the potential of the action (F.2) evaluated on 



the surface ?/> r = 0, and the target space metric is just the metric induced on the surface from 
the metric qAB- 

To obtain the higher order corrections we can proceed as follows. The equation of motion for 
the fields ip v is 



Vt/> T + 2r^ e V X a • W w = H'^ a V a + H o Mq^V - (FA) 



1,6 



2„£ai 



Here the connection coefficients are those of the target space metric qAB- We next expand this 



equation to linear order in ip r and use the condition (F.3) to obtain 



+2rf e v x a • vv 



e 



Hfa^^Va - M 2 q sr V re 

2„£<n 



(F.5) 



where all the metric coefficients, connection coefficients and their derivatives are evaluated at 
ifj T = 0. Now in the large M or adiabatic limit, the dominant term on the left hand side will be 
the term proportional to M 2 , and dropping the other terms gives a simple algebraic equation for 
the leading order contribution to ip r : 



[q^Vre] ^ 



1 

w 



(F.6) 



Substituting the solution given by Eq. (F.6) into the action (F.2) gives the required, 0(1/ M 



corrections to the action. The first term on the right hand side of Eq. (F.6) will give nonlinear 



corrections to the kinetic energy. (We assume that the second fundamental form or extrinsic 
curvature of the surface ^ r = is nonzero, otherwise these corrections would vanish.) 
As a simple example, consider the theory 



C 



1 



1 



1 



(Vxr - ; W - o M > + V<(Vx)7"v 



(F.7) 



The equation of motion for ip is Dtp — M 2 ip = (Vx) 2 /m p with leading order solution ip 
— (X7 x) 2 / ( m pM 2 ) . The corresponding corrections to the action for x scale as (Vx) 4 /(m-pM 2 ), 



m 



agreement with Eq. (5.6). The scaling (5.6) of other operators can be derived similarly. 
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